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We consider warped compactifications in (4 + d)-dimensional theories, with four dimensional de 
Sitter dSi vacua (with Hubble parameter H) and with a compact internal space. After introducing 
a gauge-invariant formalism for the generic metric perturbations of these backgrounds, we focus on 
modes which are scalar with respect to dSi. The physical eigenmasses of these modes acquire a 
large universal tachyonic contribution — l2d/(d+2)H 2 , independently of the stabilization mechanism 
for the compact space, in addition to the usual KK masses, which instead encode the effects of 
the stabilization. General arguments, as well as specific examples, lead us to conjecture that, for 
sufficiently large dS curvature, the compactified geometry becomes gravitationally unstable due to 
the tachyonic growth of the scalar perturbations. This mean that for any stabilization mechanism 
the curvature of the dS geometry cannot exceed some critical value. We relate this effect to the 
anisotropy of the bulk geometry and suggest the end points of the instability. Of relevance for 
inflationary cosmology, the perturbations of the bulk metric inevitably induce a new modulus field, 
which describes the conformal fluctuations of the 4 dimensional metric. If this mode is light during 
inflation, the induced conformal fluctuations will be amplified with a scale free spectrum and with 
an amplitude which is disentangled from the standard result of slow-roll inflation. The conformal 4d 
metric fluctuations give rise to a very generic realization of the mechanism of modulated cosmological 
fluctuations, related to spatial variation of couplings during (p)reheating after inflation. 



I. INTRODUCTION 

De Sitter or quasi-de Sitter 4d geometries describe the 
present day acceleration of the Universe as well as the 
inflationary expansion at very early times. This stim- 
ulates significant interest towards the construction of 4 
dimensional dS geometry in the context of fundamen- 
tal string/M theory, which are formulated in 10/11 di- 
mension, with a compact internal space Q, 0, Q ■ Sim- 
ilarly, compactification to 4d dS geometry takes place 
in phenomenological braneworld models, where the inner 
space is periodic with orbifold branes at the edges 0. 
Although most of the activity in this area has been en- 
hanced by very recent progress (both on the observational 
and theoretical side) the issue of dimensional reduction 
to the outer cosmological space-time was popular since 
the 1980s, either in high dimensional supergravity theo- 
ries or on phenomenological grounds. For example, see 
the collection of references on KK cosmology given in [f| . 

The geometries we discuss here have d spatial dimen- 
sions wrapped on a compact manifold M. , in addition to 
the standard (3 + 1) space-time. Many of the mentioned 
examples are covered by the (4+G0 dimensional geometry 
with the metric 

dg 2 = e 2A(y) (_ M 2 + e 2Ht d ^ + e ~^V) fa ( y ) dy a dy b ^ 

where the outer space is a dS 4d metric with a Hub- 
ble parameter H, while g ab is the metric of the compact 
inner space with d coordinates y a . For generality, we 
also include the warp factor A(y). In the following we 
use greek letters to describe the (3 + 1) outer space-time 
coordinates (fi = 0, ...,3) while roman letters span the 
inner compact space coordinates only. Capitalized ro- 



man letters span all coordinates. With the re-definition 
of the inner space metric g a b = e~ iA ( y ' g a b and the warp 
factor A(y) = e A ^ v \ we can re- write in a conformally- 
factorized form 

ds 2 =A{yf{dsl+g ab dy a dy b ) , (2) 

where ds\ is de Sitter metric. This form of the metric is 
more convenient for developing the formalism of metric 
perturbations around a warped dS& x M. background, 
which we will present below. 

^From the string theory/supergravity perspective, re- 
cent studies have concentrated on the compatibility of 
high (4 + d) dimensional geometries (d = 6) with a 4d 
de Sitter geometry as the outer space-time, and on the 
stabilization of the internal space to dS± • The bulk ge- 
ometry, which is usually treated in the supergravity limit, 
requires a careful study of the (4 + d) dimensional (bulk) 
Einstein equations. The progress in finding solutions is 
related to the identification of various possible sources for 
the bulk stress energy tensor Tg , including supergravity 
lagrangian fields, branes with fluxes, etc. In phenomeno- 
logical braneworld models, the (4 + 1) or (4 + 2) dimen- 
sional bulk geometry is stabilized, for instance, by means 
of bulk scalar fields as in the Goldberger- Wise model [|| , 
or more generally through the Casimir effect 0, II] • 

Very often, however, the stabilization is studied at the 
level of a 4 dimensional effective theory, where the in- 
ner space geometry emerges in terms of moduli fields. If 
the effective potential of the moduli has a minimum, the 
moduli are considered to be stabilized. An example will 
be the stabilization mechanism with quantum field the- 
ory effects which emerge from the properties of M. , e.g. 
the Casimir effect from compact inner dimensions. It is a 
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reasonable assumption that in principle there shall be a 
bulk prototype (Tg) (may be as a complicated functional 
of the metric and the topology) which generates corre- 
sponding terms of the 4d effective potential. The version 
of string theory dimensional reduction of [j| is realized 
with the participation of the instanton effects (which are 
manifested in 4dim superpotential of N = 1 supergrav- 
ity). An interesting question is whether it is possible to 
think about a bulk prototype Tg of these effects. 

While it is reasonable to assume that the low energy 
4d effective action is Einstein gravity plus moduli fields 
(or Brans-Dike gravity plus moduli), this picture lacks 
the full higher dimensional solution of the Einstein equa- 
tions with proper sources. In particular, the four di- 
mensional description is inadequate for studying the high 
energy regime (high dS curvature) where the stabiliza- 
tion can break down. For example, in the 5 dimensional 
braneworld models, the issue of stabilization due to the 
bulk scalar field can be formulated fully in terms of the 
stability of the 5 dimensional warped geometry against 
scalar metric perturbations. From this study one can de- 
termine the modes which enter in the exact low energy 4 
dimensional description, namely the radion and the other 
KK modes of the system |9| , which cannot be found from 
an heuristic 4 dimensional effective potential. 

We are not concerned here with the details of the sta- 
bilization mechanism. Our main goal is to determine the 
effects of the dS expansion on the stability of these ge- 
ometries. As we will show, compactifications to a 4dim 
dS space-time are more difficult to achieve than com- 
pactification to a 4dim flat space-time. This is due to a 
tachyonic contribution to the square of the mass of scalar 
metric perturbations arising from the de Sitter curvature. 
In what follows we will discuss the properties of (4 + d) 
dimensional classical Einstein equations assuming a bulk 
Tg as a source, but without specifying it. The gravita- 
tional instability which we will discuss comes from the 
gravity sector, so that the exact form of Tg will not be 
crucial. Of course it will be attractive to check how the 
effect works for each particular model of Tg. However, 
we will try to argue that the instability effect is generic. 

Previous studies of the gravitational sector (i.e. met- 
ric perturbations) of the compactification concentrated 
on the dimensional reduction to 4dim flat space-time p| . 
The instability effect which we find emerges when the 
outer space is curved instead of flat, and it is propor- 
tional to the curvature. Notice that test scalar fields 
propagating in dS& x M. do not exhibit tachyonic masses 
under KK projection Thus, there is a significant dif- 
ference in the mass spectrum between the test fields and 



the self-consistent treatment of the metric fluctuations 
around this background. 

The plan of this paper is as follows. In Section [n] 
we introduce a generic formalism of metric fluctuations 
around the dSi x M. background and define scalar, vec- 
tor and tensor fluctuations acording to their transforma- 
tion properties relative to dS± . In Section II I II we intro- 
duce the linearized equations that govern the evolution 
of scalar perturbations of the metric and show how the 
tachyonic instability of the scalar modes arises from the 
gravitational sector. In Sections IIVI and 0we show two 
explicit example where the instability is manifested. In 
Section lyTl we discuss how the instability may be relevant 
for the generation of the primordial perturbations, pro- 
viding a very natural and model-independent realization 
of the idea of modulated perturbations |25l |26| . Conclu- 
sions and a brief discussion can be found in Section [VIII 



II. METRIC PERTURBATIONS AROUND THE 
BULK GEOMETRY 

The setting of the problem we are discussing is simi- 
lar to the routine setup for (3 + 1) dimensional cosmol- 
ogy with the metric ds 2 = a(t) 2 (—dt 2 + da 2 ), where da 2 
is the metric of 3d space of constant curvature. Per- 
turbations of this metric can be classified in terms of 
freely propagating gravitational waves, together with vec- 
tor and scalar perturbations, supported by the perturba- 
tions of the sources (for instance scalar field fluctuations 
5(j) in the scalar field cosmology). The classification of 
cosmological fluctuations is based on the gauge-invariant 
Bardeen's variables 0] (see also [l3^- 

In the background geometry there are also different 
modes of metric perturbations, corresponding to differ- 
ent irreducible representations of the background symme- 
tries. Modes living in different representations are not 
coupled at the linear level. For this reason, it is very 
useful to extend the classification of 0, to this case 
too. 



A. General formalism of (4 + d) dimensional metric 
fluctuations 

Bearing in mind the similarities in the structures of 
our bulk metric J5J) and the (3 + 1) cosmological metric, 
the generalized perturbed form of the metric can be 
written as 



ds 2 = A 2 (y) {g ac (S c b + 2cj>$)dy a dy b + [(1 + 2^)g^ + 2E tliV + 2F {ll . v) + V] ^dx v - 2A a ^dy a dx»} (3) 

I 

with g liv dx^dx v = ds 2 and where and h^ v are di- A afl can be further decomposed into individual scalar and 
vergence free and /i M „ is traceless. The cross component 



3 



vector contributions as A all = B a -i + S a i- It is simple to 
check that the variables r/>, E, F^, A afll account for 
all possible physical metric perturbations. 

Gauge invariant variables that are counterparts to the 
Bardeen variables in (3 + 1) dimensions can be derived in 
this setup by considering infinitesimal coordinate trans- 
formations x A — ► x A — x A + ^ A . In particular we have 
new scalar variables 

$ab = 4>ab + ^-{B c -E, c )g ab -{B (a -E { . a ). b) ,{A) 

* = i/>+^-(B c -E. e ). (5) 

In terms of jour dimensional tensor properties in the 
4d observable world, we distinguish scalar perturbations 
4> bl ip, E, vector perturbations F lll A all and tensor pertur- 
bations hfj, u . 

This is a rather important and subtle point. What is a 
scalar fluctuation with respect to 4 dimensional de Sitter 
space-time, may have a non-scalar nature in the bulk 
geometry. For instance, a free graviton in the bulk may 
have a gravi-scalar projection in the 4d space-time. Be- 
low we will construct a specific example where the scalar 
metric perturbation of Q in fact comes from free bulk 
graviton. 

Our grouping of the metric perturbations J2J is very 
similar to the metric fluctuations formalism of | 14| de- 
signed for higher dimensional braneworlds. The impor- 
tant distinction between our formalism and that of 0| 
are the following: we keep the time coordinate in the 
homogeneous part of the de Sitter metric, and the tenso- 
rial properties of fluctuations are related to the 4 dimen- 
sional symmetric part. In 14] the tensorial properties 
are related to the 3 dimensional homogeneous space. In 
addition; we do not rely in any way on the existence of 
branes in the models. On the other hand, our formal- 
ism can embrace the braneworld inflation as a particular 
case. 

As mentioned, the main advantage of using the split 
© is that perturbations of different types are not cou- 
pled at the linear level. The analysis is further simplified 
if the same decomposition is performed also for the per- 
turbations of the source term 8Tg . In this work we study 
metric perturbations of the scalar type. [30l | 

B. Scalar perturbations 

For this case we can define a generalized longitudinal 
gauge analogous to the (3 + 1) dimensional case. In this 
gauge the scalars B a and E vanish, and the scalar per- 
turbations <pab and if) coincide with the gauge invariant 
variables Q a b and 'J, respectively. Thus, without lost of 
generality, the 'scalar' part of metric fluctuations around 
a geometry with dS± x M. symmetry can be simply writ- 
ten as 

ds 2 = A 2 (y) {g ac (5l + 2^ b )dy a dy b + (1 + 2^4} (6) 



We can obtain a further reduction in the degrees of 
freedom where the internal space M is maximally sym- 
metric. In this case (f> b — * 5 b <fi i and we obtain 

ds 2 = A 2 (y) [(l + 2^(x,y))ds 2 + 

+ (l + 2<S>(x,y))g ab (y)dy a dy b ] . (7) 

Eq. Q is the the starting point for our treatment of 
scalar perturbations It is instructive to compare the sim- 
ple form of scalar metric fluctuations in the longitudi- 
nal gauge around a (4 + d) symmetric dSi x M. back- 
ground with lower dimensional examples. If we replace 
A(y) -> a(t), ds\ -> da 2 , g a bdy a dy b -> —dt 2 , eq. re- 
duces to the familiar form of scalar metric perturbations 
around a (3 + 1) dimensional FRW cosmology. Secondly, 
if g a bdy a dy b = dy 2 is simply one extra dimensional line 
element, eq. Q gives us scalar metric perturbations in 
the (4 + 1) braneworlds 0,0, where $ is related to the 
radion. 



III. LINEARIZED EQUATIONS FOR SCALAR 
PERTURBATIONS 

First we perform a customary perturbative analysis of 
the linearized Einstein equations for the scalar metric 
perturbation Q. It is remarkable that, for symmetric 
A4, the two scalars $ and \f completely describe all the 
scalar metric fluctuations around dS x M. . The modes $ 
and ^> will be coupled to scalar perturbations of 5T- , if 
they are present. Typical example of scalar fluctuations 
in 5Tg are the fluctuations of the bulk scalar field, which 
is used in the stabilization of M. . 

We will show that $ and ^ are proportional to each 
other, and decompose ^ into scalar harmonics (KK 
modes) on dS^. We will then construct the general ex- 
pression for the lowest eigenvalue (square of the physical 
mass) of the scalar harmonics, and show that it contains 
a tachyonic contribution, which signals a possible gravi- 
tational instability. 

It is very important to stress that <I> and ^ exist as dy- 
namical fields even if STg is absent, contrary to the case 
of (3 + 1) dimensional cosmology. Indeed the bulk geome- 
try dS x M. is anisotropic as a whole. Therefore we may 
expect the presence of an anisotropic mode of the gravita- 
tional instability, even when the fluctuation of the source 
STg are absent. It turns out that this pure gravitational 
mode is also manifested in the scalar fluctuations $ and 
"J. In the next section we put the metric (7J into the 
context of the generalized anisotropic Kantowsky-Sacks 
geometries and show that the longitudinal metric fluctu- 
ations around correspond to gravitationally unstable 
modes of the anisotropic metric QJ. 

A. Linearized Einstein Equations 

Consider the linearized Einstein equations for the 
scalar perturbations <& and 'J. The equations for 5G£ 
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contain a diagonal part, proportional to 5£, plus a sec- 
ond term which is nontrivial also for ji ^ v . It amounts 
to 



SGI = -A- 2 V M V„ [2* + d$] , (i^v 



(8) 



where V^ t is the covariant derivative with respect to 
the non-compact de Sitter space. In the absence of an 
anisotropic stress tensor in the outer space, ST^ = for 
Ijl ^ v , (which is typical in cosmology) , eq. (JSJ) imposes 



4> 



(9) 



This formula replaces the well known relation for $ and 
\E' scalars in (3 + 1) dimensional cosmology, $ = — ^ [lflj . 
and generalizes the result $ = — 2 "J valid for the (4+1) 
braneworlds [l7|. 



The linearized Einstein equations for the scalar pertur- 
bations \& yield 



<5G£ = A- 2 6t: 
5G» = 



6H 2 



d + 2 



d+2 



R + L l 



<5G£ 



d+2 



r 



* = 8ttGST^ oc <5£ , 



2D, 



is 



4- = 8ttG (5T,f 



(10) 

(11) 
(12) 



In these equations D 4 = —d 2 — 3Hd t + e~ 2Ht d 2 is the 
scalar V M V M operator of the (3+1) dimensional de- 
Sitter space, while V a , G£ = i?^ — ^5% R, are com- 
puted on the compact manifold M. with the metric g a t ■ 
LxjL^Ls denote terms proportional to derivatives of A . 
Such terms do not contain second derivatives with respect 
to a+ and can be neglected in the present discussion. For 
completeness, they are reported in Appendix. Finally, G 
denotes the higher dimensional Newton constant. 



B. Tachyonic Contribution to the Eigenmasses 

The scalar perturbations as well as the scalar per- 
turbations of the stress-energy tensor 5T- , due to the 
symmetries of the background geometry can be decom- 
posed into the scalar eigenmodes Q n (x) 

9(x, y) = 5^* n (i0 Qn(x) , 

n 

6T£(x,y) = Y^^BrM Qn(x) , (13) 

n 

where the sum corresponds to the the KK tower of the 
4 dimensional scalar modes. These modes Q n , which are 
common for all fluctuations, obey the free massive scalar 
field equations 



-□ 4 + m 2 n ) Q n = 0, 



(14) 



where m 2 are the separation eigenvalues of the decompo- 
sition jn§. Substituting (JT3J and (O in (T2), we obtain 
a set of equations for each eigenmode. In practice, we 



simply replace 5T£) by {# n , £Z# n } , and D 4 by 
the corresponding eigenvalue m 2 . 

Among l|10|) - (|12|) . l ne equation (|12|l is the dynamical 
equation while equations UlUfl . (|11|) are the constraint 
ons. Moreover, these equations are connected by the 
Bianchi identity. Thus, it will be enough to work with 
the two equations i|ll[) and (|12fl only. 

Let us consider the dynamical equation (|12fl . We make 
the following observations which is the crucial one for 
our discussion. We note that the four dimensional op- 
erator only enters in the diagonal components of the 
Einstein tensor 6G% corresponding to the internal space, 
see eq. (|12|l . Most importantly, it always appears in the 
combination EI4 + -j^H 2 - Hence, this combination will 
always generate the contribution 



□4 



\2H 2 d 



= m n + 



2 

12 H 2 

1 + 2/d 



u 2 \P 2 



(15) 



in the equations for the modes along M. . This is true 
for any choice of M. and for any underlying stabiliza- 
tion mechanism, since it is due to the symmetry of the 
background configuration Thus, the eigenvalue m\ 

from the outer space-time always appears together with 
12 H 2 

1+2/d' 

After the factorization l|13|l . the resulting equations are 

the equations for the functions j^nd/), ^Tg n {y)\ on 

M . Typically, they can be reduced to a second order el- 
liptic equation. For compact Ai, this corresponds to the 
eigenvalue problem which determines the combinations 
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n\ . Once these values have been obtained, the physical 
masses m 2 n of the perturbations immediately follow from 
eq. lfT5| . 

Further we can use the symmetry of M . If the warping 
is absent, A = const, on general ground we expect that 
the solution of the eigenvalue problem on A4 will depend 
on the curvature of A4, [i^ ~ R; we denote it as ix n (R) 2 . 
On the other hand, for non-vanishing warping operators 
L may bring in ^ extra dependence from A(y) in terms 
of its characteristic scale, which we denote as fi n (A) 2 . 

Taking all this together, the eigenvalues of the scalar 
metric perturbations in dS space-time are expected to 
have the following structure 



i , jp- . The late-time asymptote of this solution 



12 H 2 

T+2/d 



fj? n (R) + » 2 n (A) 



(16) 



To support this conjecture, in the next Sections below 
we will give examples of explicit computations of m\. 

The most important conclusion is that the curvature 
of non-compact de Sitter dimensions gives a negative 
(tachyonic) contribution in eq. (|16fl to the eigenmasses 
m 2 n of the scalar fluctuations. If de Sitter curvature ex- 
ceeds a critical value so that becomes negative, we 
encounter gravitational instability of the dS± x M. geom- 
etry. This result generalizes a similar finding obtained 
in [l7| in the case of (4 + 1) dim braneworlds, where 
= -AH 2 + (J-n{A). In [lTj the interpretation of the 
tachyonic instability in m\ was not given. Here we ar- 
gue that this is an effect of gravitational instability of 
anisotropic geometries dS^ x M.. The present discus- 
sion shows that the instability is a generic issue that 
is not confined solely to the braneworld models studied 
in [TtI ITsj . Rather, it is a general property of the geome- 
tries of the form Q , which includes a very wide class of 
theories with extra dimensions. In the next Section we 
discuss how this instability emerges in different contexts. 



C. Tachyonic Instability of dS± x M 

We now follow what happens if m 2 n is negative. For this 
we have to turn to equation (|14|) for the four-dimensional 
eigenfunction Q n (t,x). The four-dimensional massive 
scalar harmonics Q n can be further decomposed as 
Q n (t,x) = J f { h n) (t)e^ s d 3 k. The temporal mode func- 
tions /£ (t) obey the equation 



f + 3Hf+ (e- 2Ht k 2 + ml) f = 0, 



(17) 



where dot denotes time derivative, and we have dropped 
the labels k and n for brevity. We choose the so- 
lution of (|17fl which corresponds to the positive fre- 
quency vacuum-like |3l| initial conditions in the far past 
t -> -oo, f k (t) - ^e^, r) = fdte~ m . For the 

tachyonic mode (m^ < 0), the solution to equation l|17|) 
with this initial condition is given in terms of Hankel 
functions f ( k n) (n) = ^-H^l 2 ^ (krj), with the index 



' 9 | \ m n I 

~~H T ~ 

diverges exponentially as t — > oo (n — > 0) 
with the numerical factor 



A,, 



'9 Kl 3 
A H 2 2 



(18) 



(19) 



The leading term is associated to the lowest negative 
7Xq ~ —few x H 2 , which defines the 



mass square m n — 
instability of metric fluctuations 

^(t,x;y) oc e XoHt 



(20) 



In the following sections we consider examples where 
ml can be calculated explicitly, and we will provide an 
interpretations of the tachyonic instability. 



IV. EXAMPLE: INSTABILITY OF (4 + 1) 
INFLATING BRANEWORLDS 



The inclusion of additional fields is typically required 
to achieve the stabilization of the internal space. How- 
ever, this complicates significantly the study of the per- 
turbations of the system, so that explicit results have 
been obtained only for the simplest configurations. The 
theory of metric perturbations was developed for (4+1) 
Braneworld models with inflation. Although we do not 
deal specifically with the braneworlds in this paper these 
provide a good example for models that are affected by 
the results discussed in this work and constitute an in- 
teresting set-up where calculations have been performed 
in detail. Here we consider the (4 + 1) dimensional 
braneworld with de Sitter branes of codimension one, 
with and without the presence of a bulk scalar field. The 
geometry of the system, with the inclusion of scalar per- 
turbations, is again of the form (JJJ, with parallel branes 
as edges of the internal coordinate y . The function A (y) 
is known as warp factor. 



A. Inflating Braneworld without stabilization 

First let us consider an inflating braneworld without 
stabilization (i.e. without bulk scalar field). We assume 
a negative cosmological constant in the AdS bulk. Scalar 
metric perturbations are given by the form (J7J, and from 
@ we have $ = —2^. The Einstein equation for the 
linearized perturbations $ has to be supplemented by 
the junction conditions at the orbifold branes, for details 
see H, [13,113 ■ Equation JTIJ gives us 



A 1 

*' + 2— * = 
A 



(21) 



where prime stands for d y . This relation shows that 
the KK tower is absent in this case, and that the only 



6 



eigenmode ~ I /A 2 is present. Equation (|12(l gives 
(□ 4 + AH 2 ) ^ — 0. We then see that he perturbation ^ 
is described by the free wave equation. In other words, 
for the branes without stabilization, \P is nothing but the 
gravi-scalar 4 dim projection of the bulk gravitational 
wave. Second, the square of the mass for this mode is 
tachyonic 



= -AH 2 



(22) 



in exact agreement with our eq. I(I6() . In this case 
Hn{A) = 0, (J>n(R) = 0. The braneworld between two 
de Sitter branes is unstable. The equation 1)22(1 for in- 
flating branes without stabilization was obtained in |l9j| . 



B. Inflating Braneworld with Bulk Scalar 
Stabilization 



The addition of a bulk scalar field can, in principle, 
provide the stabilization of the braneworld configuration. 
The scalar field <f> can have a potential V in the bulk and 
potentials U% t 2 at the two branes. The internal coordi- 
nate can be stabilized by suitable choices of these poten- 
tials. However, as we will see, the branes can be curved 
to such a degree that the stabilization mechanism will 
fail. 

Let us consider a stationary configuration, character- 
ized by a fixed interbrane distance, and let us study its 
stability. Perturbations of the bulk field 5<fr are cou- 
pled with the scalar perturbation of the metric, as we 
described above. We can decompose them as in eq. (|T3)l , 



' (x, y) = Yl ^« ( y ) ® n ( x ) 



(23) 



Next, we can use the well known expression for the bulk 
stress energy tensor 8Tp in terms of 8<f> and the back- 
ground field 4>(y) (see [H], for details), and substi- 
tute them in the linearized Einstein equations. The con- 
straint equation l|l(J|) gives a connection between ^ n (y) 
and 5<j) n (y) 



(A 2 * n y 



A 2 



Equation (|12|l yields instead 



4? (A 
3A 



(24) 



(25) 



Here, combinations of A and its derivatives are expressed 
in a more compact manner through </>' ', using the back- 
ground equations. 

Because ^ and 6(f) are related by the constraint equa- 
tion (|24l) . in fact there is only one dynamical degree of 
freedom, which is linear combination of the pair. The 
mass of the four dimensional modes is most easily ob- 
tained from an eigenvalue equation for the wave functions 
of the modes along y . 



Combining equations 1(24(1 and 1(25(1 one obtains 

dj ry 1-) £ 2 

In ^ (A 3/2 



AH 1 - 



(A 4') 



Un = 



-1/2 



(26) 



The mass spectrum of the system is determined by the 
eigenvalue problem, i.e. by the solutions of this equation 
which satisfies the boundary conditions for the perturba- 
tions at the two branes. It can be shown that the lowest 
eigenmass is ^3 



-AH 2 + fM) {A) , 



(27) 



where (i (A) » 3 //J"^ 
stein equations we have </ 



From the background Ein- 



G 



A' 2 



,A" 



3H 2 



that (J.q(A) is expressed through the warp factor and its 
derivatives. 

As we have argued in the general case (set d = 1 in 
eq. 1(15(1 ). the eigenmasses m n enter in the Einstein equa- 
tions in the combination m} n + A H 2 . Hence, the 4 H 2 
term results in a negative contribution to the physical 
masses of the four dimensional modes. The stabilization 
mechanism has to be "sufficiently strong" to counterbal- 
ance this effect, and the configurations for which the right 
hand side of eq. ((27(1 is negative are subject to tachyonic 
instability (|2Til) . 

Since the warp factor and the expansion rate H both 
enter in /j,q (A) , this term should be considered as a func- 
tional of H . Hence H affects m 2 , through both terms. 
However, no (A) only mildly depends on H and by far, the 
strongest H dependence is given by the negative —4 H 2 
contribution. Therefore the critical value of dS curvature 
H = H c which borders stable and unstable brane configu- 
rations can be estimated from AH 2 = /zo(A), where ^o(^) 
can be calculated explicitly for concrete braneworld mod- 
els (for specific examples of the calculation of no(A), see 

MM)- 

The study of the linear perturbations around the sta- 
tionary configuration gives information about its stabil- 
ity, but it cannot be used to determine the whole dy- 
namics of the system when the configuration is unsta- 
ble. The dynamical evolution can be studied numerically 
with the branecode package, a numerical code [Tsl ] de- 
signed specifically for this task. It was found numeri- 
cally (by methods which are completely different from 
the linear analysis described above) that indeed station- 
ary inflating branes configurations with too high H are 
unstable. Moreover, the full numerical treatment of the 
dynamics reveals the end points of the non-linear evo- 
lution. It turns out that, depending on the parame- 
ters, unstable brane configurations can be re-structured 
to stable inflating branes configurations with lower H 
(H < H c ).\i the parameters of the branemodel do not 
admit another stationary configurations, the branes are 
instead colliding. Bulk geometry of the collapsing branes 
asymptotically reaches the 5 dim Kasner solution ds 2 = 
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~dt 2 + t^dy 2 + ^Li t 2pi dx 2 . Three of the Kasner in- 
dexes Vi are the same, as disctated by the brane isome- 
tries fl8l. 



V. EXAMPLE: INSTABILITY OF dSi x S d 

As the next example we consider a (4 + d)-dimensional 
space-time with the internal dimensions compactified on 
a d— sphere of radius r$ . This is a toy model and has 
no application to compact extra dimensions. However 
it serves as a simple geometry with which we can check 
our conjecture from the GR point of view. We work out 
in details the case in which the only source term is a 
cosmological constant A . The system is rather simple 
and its dynamics can be studied exhaustively (for the 
low dimension case of (IS2 X 62 see He}): there is one 
stationary configuration of the type Q), which, due to 
the lack of a stabilization mechanism, turns out to be 
unstable. Our main point is that the instability can be 
computed precisely with the general approach outlined in 
Section IlIII Below we will apply this generic formalism 
to the case of dS& x Sd with a positive bulk A . 

However, there is a complementary treatment of sta- 
bility issue of dS/t x Sd based on the investigation of 
its gravitational evolution within a more general class 
of anisotropic geometries, the Kantowsky-Sacks metrics, 
of which dS*4 x Sd is a particular solution. We show 
that the two approaches are equivalent. This provides 
insight on the nature of the tachyonic instability conjec- 
tured for dSi x M. , as the general gravitational instability 
in the class of anisotropic metrics. The dynamics of the 
the Kantowsky-Sacks metrics in (4 + <i)-dimcnsions can 
be analyzed by the qualitative methods of the dynamical 
systems The phase diagram can be extended far from the 
stationary configuration, so we can find the end points of 
the tachyonic instability. 



A. Instability of linear scalar perturbations around 

dSi x Sd 

We consider the simple situation in which the geometry 
is factorized (A = 1 in eq. J7J) and the internal manifold 
M is a sphere of dimension d and radius r . We assume 
that the only source term is a cosmological constant A. 
This systems admits a stationary solution with static Sd 
and de Sitter external space. This configuration is char- 
acterized by 



R = MH 2 



H = 



2A 



3(d + 2) 



(28) 



The curvature of the inner space is R = d(d — l)/^. 
The system allows for expending H > and contracting 
H < phases. 

Consider linearized scalar perturbations around 
this background in the absence of any dynamical field in 



the system ST^ = 0. The first Einstein equation l|10|) is 
satisfied identically due to relation l|28|) . The constraint 
equation (fTT|) gives us \P = const, so that only the light- 
est mode is present; while the tower of KK modes is ab- 
sent (this is similar to the case of a braneworld with neg- 
ative bulk cosmological constant without stabilization; 
also in this case the KK tower would be recovered upon 
introduction of dynamical fields). The dynamical Ein- 
stein equation (|12|l gives us 



□4 



12d 
d + 2 



H 2 - 



4(d-l) 
d(d + 2) J 



R * = 



^From here, we immediately have 



12 H 2 

1 + 2/d 



(29) 



(30) 



where Hq{R) = fj^+fy-R- The term ^{R) , proportional 
to the curvature of the d— sphere, has comparable mag- 
nitude and opposite sign with respect to the tachyonic 
term. However, the net mass squared of "J is negative for 
any dimension d, 



= -6H 



(31) 



This signals the tachyonic instability of the configura- 
tion l|28(l . Interestingly, the ratio mo /H does not depend 
on d. 

Applying eq. (|20[1 to this example, the linear scalar 
perturbations are growing with time according to 



^{t,x;y) oc e 



\ Ht 



exp 



33-3 



Ht 



(32) 



Notice that l)29|) is the free wave equation; this means 
that corresponds to the gravi-scalar projection of the 
bulk gravitational wave. 



B. Gravitational instability of anisotropic 
geometry dS& x SV. 

The expression Ij32(l controls the development of the 
gravitational instability as long as l^l <g; 1 . However, to 
understand the nature of the instability, we have to study 
the dynamical evolution beyond the linear regime. To do 
so, we allow the radius of M. to be a time dependent 
function r (t) , and we consider a general FRW evolution 
for the non-compact space, 



ds 2 = -dt 2 +a 2 (t) d 3 x+r 2 (i) [d$l + sin 2 6»i (d6 2 + ■■■)] , 

(33) 

where 9%, 0%, . . . are angular coordinates on the d— sphere. 
For definiteness, we specify to d — 7 (the evolution of 
the system in a lower dimensional context was studied in 
details e.g. in |2(|). From the (0,0) Einstein equation, 
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the expansion rate of the non-compact dimensions can 
be written as a function of r and its time derivative r , 



a _ -21r ± [4Ar 2 + 63 r 2 - 84] 
a 6 r 



1/2 



(34) 



We can substitute back this relation into the remaining 
equations, and obtain a second order differential equa- 
tion for r alone. As before, we can rescale A away by 
redefining x = VAr , and r = \/At . We then obtain the 
following dynamical system 



dx 

dx 
d7 



(35) 



Ax 2 + 81y 2 - 108 T $V3y [Ax 2 + 63y 2 - 
18a: 



,1/2 



whose evolution is summarized in the phase portrait of 
figure n] 

The set of trajectories lies in a nontrivial manifold, 
characterized by a forbidden region (Ax 2 + 63y 2 — 84 < 0, 
bounded area of fig. QJ, and by two branches which are 
connected at the boundary of this region but discon- 
nected anywhere else. The two branches are shown in 
the right and left part of fig. ^ with the right (left) part 
corresponds to the upper (lower) signs in eqs. (|34|) . (|35|l . 
Trajectories are shown by continuous lines, while the two 
dashed lines separate the region for which the external 
space is expanding (H > , lower part in both branches) 
from the (upper) one for which it is contracting (H < ) . 
The trajectories which reach the boundary of the for- 
bidden region in the right branch are identified without 
discontinuity through opposite points along the bound- 
ary into the left branch. The background configuration 
dSi x Sj studied in the previous Subsection corresponds 
to the saddle point in the right branch, characterized 
by a static radius {x = y/27,y — 0} (in agreement with 
eq. (j2Hl); and by a de Sitter outer space. The region of 
the phase portrait around the saddle point 5+ is zoomed 
in the right panel of fig. ^ 

The linearized system l|35[) around this point is 



This is precisely the result (I32|) for "J obtained with the 
previous calculation. This confirms the equivalency of 
the two approaches for the study of the instability of the 
stationary configuration. 

The eigenvalues of the matrix l|36[> are of opposite sign, 
indicating that the stationary configuration is a saddle 
point. Its eigenvectors describe the critical directions in 
the linear regime close to the saddle points. The motion 
along the unstable line leads to two very different asymp- 
totics along the two opposite directions. Starting with a 
slightly larger r the system goes toward the decompactifi- 
cation of M. , with the contribution from the cosmological 
constant dominating more and more over the curvature 
of the d— sphere. The asymptotic regime is a straight 
line in the phase portrait, y/x = constant. Substitut- 
ing this "ansatz" in eq. Q35|). one finds y/x — > 1/V45, 
as x, y — + oo . This solution corresponds to a de Sitter 
expansion of the internal space. From eq. i|34|) we see 
that the non-compact space also approaches the same 
de Sitter expansion, 



a r 
a ' r 



(38) 



Hence, the system tends toward dS\\ , with an isotropic 
expansion of all the coordinates. Notice that the bulk 
expansion rate l|38|) is slower than the initial expansion 

rate of the outer space H = \J~^f- Isotropization of ex- 
panding solutions is produces by the bulk cosmological 
constant. 

Alternatively, the instability can lead to the collapse 
of the internal manifold. In the phase diagram of fig- 
ure n this possibility is described by trajectories which, 
starting from the stationary point, go from the right to 
the left branch, and then continue toward r — ► , while 
the external space expands with increasing expansion 
rate. The asymptotic evolution is of the Kasner type 
ds 2 = -dt 2 + t^dx 2 + t^dn 2 .. We have 



a cx (t* — t) 



V21-1 



r oc (t* —t)' 



(39) 



d 

d7 



8x 




(36) 



Next, we have to diagonalize he matrix (|36|l . The di- 
agonal components of it are diag ^ -v^-Vz^ -Vg+V^ , 

One of the eigenvalue is positive, another is negative, how 
it shall be for the saddle point. Instability corresponds to 
the positive eigenvalue. More precisely, metric instabil- 
ity growth with time as exp ( ~^®+^ 22 \/Ht ) . Recall that 



VX = H in the case of 6*7. Therefore, the metric 
instability around saddle point grows as 



where t* is the time at which the system becomes singu- 
lar (one can easily verify that l|39|l solve the dynamical 



equations at leading order as t — 
are the same in two groups, pi 
and Pi = .. - - 3^21+7 



P\\ 



70 



f* ) . Kasner exponents 

The Kasner geometry 
2ll . already advocated 



is a generic collapsing solution 

in string cosmology (23 and braneworld cosmology [Isj 
It is characteristic of a strong gravity regime, where the 
presence of any source (the cosmological constant, in the 
case at hand) can be neglected. 



VI. MODULATED PERTURBATIONS 



//V33-3, TT ' 
r(t) — Tq oc exp ( ) H t 



(37) In the previous Sections we discussed how the de Sitter 

expansion can affect the stability of the internal space. 
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FIG. 1: Phase portrait for the dynamical system 1351 . Two branches are present, according to different signs in eqs. (I34L 
1351 1. The branches are connected along the boundary of the oval region symmetrically across the axis passing through the 
point Po (as shown by the identification of points Pi and P2 in each branch. The dashed (green) curves shown the boundary 
where the Hubble parameter for the non-compact space changes sign. The stationary configuration discussed in the main text 
is the saddle point S+. Outgoing trajectories from this point either go to a decompactification of the internal space (toward 
asymptotic dSn) or they change branch and go toward Kasner asymptotic with expanding outer space and collapsing internal 
manifold. We also show a detailed view of the critical region around the saddle point S+. 



As we discussed, the stability cannot be reached by 
gravity alone. A point of phenomenological relevance 
is that the field(s) responsible for the stabilization un- 
avoidably introduce (quantum) fluctuations which per- 
turb the background geometry. Even if these fluctuations 
are heavy today, so that they have not been yet observed 
at accelerators, they could have had a mass smaller than 
the expansion parameter H during inflation. This is par- 
ticularly true if the the expansionary rate itself gives a 
negative contribution to their masses, as we discussed 
above. If this is the case, the fluctuations are ampli- 
fied to a classical level during the inflationary expansion. 
These fluctuations can arise independently on slow roll 
inflaton field. They can be generated even in absence of 
an inflaton field (if, for example, inflation is due to some 
metastable form of vacuum energy) since they are "sup- 
ported" by the bulk dynamical fields responsible for the 
stabilization of M. . 

As we will see Fluctuations of the internal manifold 
will lead to modulated perturbations. The modulated 
fluctuations were introduced recently as an alternative 
mechanism of generation of pimordial cosmological per- 



turbations |25ll26| . 

To see this, let us focus on the four dimensional part 
of the geometry. Induced metric at dSi is 



1^ = A 2 (y)(l + 2tf(a;,tf))ff /i „ 



(40) 



where g^ v is the usual four dimensional dS metric. Di- 
mensional reduction leads to appearance of the modilu 
field ^ in the four dimensional effective action. Consider 
only the lightest KK mode, after proper normalization 
we write 



S = 



(41) 

Such an action (|41() may appear if there are 1 + 3 di- 
mensional submanifolds which are point-like on M. (in 
which case 7^ is the induced metric on the subman- 
ifolds). More generally, we can also obtain it starting 
from the geometry(|7J and integrating over Ai ( in which 
case A and do not depend on y ). The first situation 
is typical of brane models, while the second is more com- 
mon of general KK theories. In both cases, we have in 
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mind a more general action than l|41|) . which will be rel- 
evant for the dynamics of the brane fields or of the KK 
modes. However, the choice (|41f) is adequate enough for 
the following discussion. 

The action (|41|) differs from a standard theory of in- 
flation due to the presenc of interaction term ^fT^. Per- 
turbations of ^ lead to rescaling of masses and couplings 
during inflation, untill the moduli is pinned down to 
its minimum where it has the the mass m\. The spatial 
fluctuations ^ of the geometry introduce a small spatial 
dependency in this redefinition, so that any physical mass 
scale and couplings will be slightly different in different 
parts of the universe. 

Inhomogeneities of the moduli field ^> are generated 
during inflation if its mass mo is ligther than the Hub- 
ble parameter. Now we recall the results of the previous 
section, that m\ contains tachyonic contribution. There- 
fore, an interesting possibilitu arises if the net value of 
m\ is smaller than H 2 . In this case during inflation we 
encounter generation of spatial fluctuations of ^. In the 
late universe long after inflation, when H significantly 
drops, we expect m\ be positive and large relatively to 
the temperature of the universe, si that ^ is stabilized. 

Small spatial variation of masses and coupingsm im- 
mediately after inflation, during (p)prehcating after in- 
flation, will result in a spatial dependency of the inflaton 
decay rate T 

AT 

-jT ~ * • (42) 

This fluctuation in the decay of the inflaton give rise to 
the adiabatic mode of the (standard) perturbations: after 
the reheating epoch the evolution of the perturbations in 
the various decay products follows the standard picture. 
The important distinction is that in this case the inflaton 
assumes a secondary role in defining the nature of the 
perturbations. 

After reheating ends and modilu $ is trapped in its 
minimum, spatial variatins of masses and couplings are 
erased. However, long wavelength metric perturbations 
generated during the inflaton deca ys are p reserved as the 
scalar cosmological perturbations \2& |2|| . 

There has been some activity in the computation of 
the spectrum of CMB anisotropies in braneworld models. 
However, many of these studies do not discuss the general 
set of metric perturbations in the bulk, but rather focus 
on the evolution of the standard four dimensional per- 
turbations with the braneworld configuration as a back- 
ground. On the contrary, the mechanism outlined here 
exploits the presence of new perturbations, related to the 
bulk dynamic, in a crucial way (notice that the conformal 
mode (|41fl is absent on the standard 4 dimensional case) . 
In addition, the one we have described is a very natural 
realization of the gen eral idea of modulated perturba- 
tions proposed in |2fiL l2fil |. since we provide a definite 
candidate (generally present in all the configurations (JTJ 
for the modulus field controlling T ) . 

There has been much focus recently on so-called 



"anomalies" emerging in the observations of LSS and 
the CMB [27j ranging from an intriguing lack of power 
on the largest scales j^, the possibility of "glitches" 
in the CMB spectrum and an indication of running of 
the spectral index of the primordial scalar perturbations 
(for a recent review, see 29] ), to detections of mild non- 
Gaussianity and indications of a finite topology for the 
3 spatial dimensions. The accuracy of future observa- 
tions is expected to increase greatly over the next decade 
and the hope is that if any of the anomalies turn out to 
be robust results they will help to constrain new physics 
modifying the dynamics in the early universe. 



VII. DISCUSSION 



A typical denominator in many extensions of the Stan- 
dard Model is the presece of extra-dimensions. In the 
simplest possibilities, the extra space is a compact and 
static manifold M . A stabilization mechanism is typi- 
cally required to explain why M. should remain static, 
while the 3 noncompact spatial dimensions are under- 
going cosmological expansion. In particular, significant 
activity has been recently made trying to reconcile this 
picture with a de Sitter (or quasi de Sitter) geometry for 
the noncompact coordinates. Indeed, observations are 
telling us that the expansion of the universe was accel- 
erating at very early times, and it is accelerating also at 
present. 

The present note is focused on the effects of the infla- 
tionary expanison on the stability of M. . We have argued 
that the dS curvature (in other words, a nonvanishing ex- 
pansion rate H) has typically the effect to destabilize the 
internal space, and that any given stabilization mecha- 
nism can be effective only up to a certain curvature. We 
have aimed to discuss this effect in the most general way 
possible, enlighting the consequences of making the only 
assumption that the geometry is of the dS x A4^ type, 
with an arbitrary compact space of d dimsensions, and 
possibly with the presence of a warp factor. Such a set- 
up is also relevant for string-theory, in the supergravity 
limit. 

To discuss the stability of the system, we have studied 
the most general set of perturbations of the dSi x A4i 
geometry. It is very convenient to classify the perturba- 
tions into irreducible representations of the dS^ symme- 
try group. The big advantage of doing so, is that modes 
belonging to different representations are not coupled at 
the linear level. The perturbations can be devided into 
scalar, vecor, or tensor modes with respect to the dS± 
isometries. After a general classification, we have focused 
on the scalar modes, since they are ususally the most 
relevant one for inflationary geometries, and since they 
encode the effects of the instability we want to discuss. 
The linearized Einstein equations for the perturbations 
show that the physical mass squares m 2 of these modes 
acquire a negative contribution due to the dS expansion, 
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(H) = - 



12 H 2 
l + d/2 



(43) 



This is true for arbitrary M. , d , and for any possible 
underlying stabilization mechanism. 

The presence of this tern is a signal for a possible 
gravitational instability of the system. Indeed, if the 
whole m 2 turns out to be tachyonic, the dS x M is 
unstable. Clearly, verifying whether (and at which H) 
the instability takes place is very model dependent issue 
which should be verified case by case (namely, for any 
given source Tg). Indeed, other contributions to m 2 
also dpend on H (although in most cases only indirectly, 
due to the fact that these terms depend on other back- 
ground quantities, and that these quantities are related 
to H through the background Einstein equations) . How- 
ever, the contribution (|43|l is generally present, and it has 
to be taken always into account. Moreover, the study of 
several different examples, as well as general arguments, 
lead us to conjecture that the term l|43|l is the dominant 
one at large H , and that no stabilization mechanism can 
be effective up to arbitrarly large dS curvature. 

In this work, we have presented explicit and exact cal- 
culations in two specific examples. The first of them 
is codimension one braneworld configurations with bulk 
scalar field(s). As shown in ^t|> m this case it is possible 
to derive an explicit upper bound m 2 < —4 H 2 + /j,q (A) 
(where A is the warp factor, see eq. I|27|) ) for the mass of 
the lightest eigenmode. This bound allows to determine 
at which H the system becomes unstable. The second 
example is when M. is a d— sphere, and when the only 
source term isw a cosmological constant. This system 
is well known to be unstable, and we have shown that 
the instability is precisely due to the tachyonic nature of 
the scalar modes which we have identified in the general 
calculation. 

There are also general arguments in support of the in- 
stability at high H . The first is related to causality. If, 
as in the standard 4 dimensional case, the Hubble length 
H" 1 has the meaning of a casual horizon, we should ex- 
pect that any stabilization mechanism cannot be at work 
when H becomes much greater than the inverse size of 
M. (since different "edges" of M would then be casu- 
ally disconnected) . Such a behavior can be found in 123 . 
where the cosmology of the Randall-Sundrum model |24| 
was studied under the assumption of radion stabilization. 
It was shown in |23| that radion stabilization cannot be 
imposed if the physical energy on the hidden brane is 
greater than ~ TeV 2 M 2 . The appearence of this inter- 
mediate scale is somewhat surprising, and indeed it was 
left unexplained in |23| . However, it can be verified that 
it is precisely at this energy density that H^ 1 becomes 



smaller than the distance between the two branes. |32| 

A second argument can be inferred from the dS± x Sd 
example. We noted that the instability has two possible 
end-points. One is characterized by a shrinking internal 
space, with asymptotic Kasner solution (a similar behav- 
ior was already noted in (0) in the case of brane colli- 
sions) . The second one is instead exact dS^+d , with all 
the compact and noncompact dimensions expanding with 
the same asymprtotic rate. This suggests that the insta- 
bility can be attributed to the tendency of an inflationary 
expansion to homogenize and isotropize the whole geom- 
etry, in this case by "dragging" also the compact space 
into a de-Sitter solution. This behavios, which is a well 
known feature of inflationary expansions, is encountered 
also in many Bianchi models. 

In the final part of the work, we have also commented 
on the possible role of the instability in the generation 
of primordial perturbations. The scalar perturbations 
appear as a conformal mode from the four dimensional 
point of view (for example, for an observer living on a 
3-brane, or after integrating out the internal space M ). 
This mode is a clear signature of the extra dimesnions, 
since it is absent in the standard 4 dimensional case. 
If, due to the tachyonic contribution, this mode light 
(0 < m < H) during the de Sitter stage, it will be am- 
plified with a scale invariant spectrum and with an am- 
plitude which is disentangled from the standard result of 
slow-roll inflation. One can think of a number of ways 
how the conformal perturbation could be responsible for 
the adiabatic mode in the later FRW evolution. For in- 
stance, consider a mixture of fields conformally and non- 
conformally coupled to the noncompact geometry. Only 
the nonconformally coupled fields will be sensitive to the 
conformal perturbations, so that we expect that the final 
mode perturbation will be proportional to the amount of 
these fields in the mixture. In the text we have instead fo- 
cused on a much simpler possibilitVjrelated to the mecha- 
nism of modulated perturbations j23,|2y]. The conformal 
factor can be usually rescaled away by a rescaling of the 
energy scales in the 4 dimensional theory (this is exactly 
what occurs in the Randall-Sundrum model |2J|). Hence, 
its fluctuations can be interpreted as fluctuations masses 
of particles and rates of processes in the 4 dimensional 
theory. In oarticular, this may be the origin of the fluc- 
tuations of the decay rate of the inflaton, which is at the 
basis of the mechanism of modulated perturbations. 
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APPENDIX A: OPERATORS IN THE 
LINEARIZED EQUATIONS FOR THE 
PERTURBATIONS 



The linearized equations (|12fl for \t have been given 
in the main text in terms of the three operators £1,2, 3 • 
These operators are 



In = 
U = 

U = 



IpM- 



v a a r~\A 

(d + 4)— V Q +4 — +2(d-l)- [2 

2^V 
A 



„VM- „V b A~ VMVfeA V a V fc A 
-2 — — V 6 -2— ^-V Q + 8 ^ 4- 



.4 



.4 



A 2 



A 



2(d + 3)V c ^V c + 4 



DA (^ A ) 



(Al) 



where quantities with tilde are computed from the back- background metric on y 
ground metric on x , while quantity with hat from the 



